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Ultracold atom research presents many avenues to study problems at the fore- 
front of physics. Due to their unprecedented controllability, these systems are 
ideally suited to explore new exotic states of matter, which is one of the key driv- 
ing elements of the condensed matter research. One such topic of considerable 
importance is topological insulators, materials that are insulating in the interior 
but conduct along the edges. Quantum Hall and its close cousin Quantum Spin 
Hall states belong to the family of these exotic states and are the subject of this 
chapter. 
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Ultracold Atoms, synthetic gauge fields, and topolog- 
ical states 

Quantum engineering, i.e., preparation, control and detection of quantum sys- 
tems, at micro- to nano-Kelvin temperatures has turned ultracold atoms into a 
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versatile tool for discovering new phenomena and exploring new horizons in di- 
verse branches of physics. Following the seminal experimental observation of 
Bose Einstein condensation, simulating many-body quantum phenomena using 
cold atoms, which in a limited sense realizes Feynman's dream of quantum sim- 
ulation, has emerged as an active frontier at the cross section of AMO and con- 
densed matter physics. There is hope and excitement that these highly tailored 
and well controlled systems may solve some of the mysteries regarding strongly 
correlated electrons, and pave ways for the development of quantum computers. 

The subject of this chapter is motivated by recent breakthroughs in creating 
"artificial gauge fields" that couple to neural atoms in an analogous way that elec- 
tromagnetic fields couple to charged particles |TJ. This opens up the possibility of 
studying quantum Hall (QH) states and its close cousin quantum spin Hall (QSH) 
states, which belong to the family of exotic states of matter known as Topological 
Insulators. 

QH and QSH effect usually requires magnetic field and spin-orbit coupling, 
respectively. In contrast to electrons in solids, cold atoms are charge neutral and 
interact with external electromagnetic field primarily via atom- light coupling. Ar- 
tificial orbital magnetic field can be achieved by rotating the gas (9J. In the ro- 
tating frame, the dynamics of atoms becomes equivalent to charged particles in a 
magnetic field which is proportional to the rotation frequency. To enter the quan- 
tum Hall regime using this approach requires exceedingly fast rotation, a task up 
to now remains a technical challenge. In this regard, another approach to en- 
gineer artificial magnetic field offers advantages. This involves engineering the 
(Berry's) geometric phase of atoms by properly arranging inhomogeneous laser 
fields that modify the eigenstates (the dressed states) via atom-light coupling [JTOj . 
Such geometric phase can mimic for example the Aharonov-Bohm phase of or- 
bital magnetic field. Because the dressed states can be manipulated on scales of 
the wavelength of light, the artificial gauge fields can be made very strong. For 
example, the phase acquired by going around the lattice unit cell can approach the 
order of it. This gives hope to enter the quantum Hall regime. The strength of this 
approach also lies in its versatility. Artificial gauge fields in a multitude of forms, 
either Abelian (including both electric and magnetic field) or non-Abelian (includ- 
ing effective spin-orbit coupling), can be introduced for continuum gases or atoms 
in optical lattices [[10). Recent NIST experiments have successfully demonstrated 
artificial magnetic field |[TJ as well as spin-orbit coupling [2] for bosonic Rb atoms 
and experiments are underway to reach QH and QSH regimes. 
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Figure 1: Figure illustrates distinction between ordinary band insulator (left), 
QH (middle) and QSH (right) insulators in the energy spectrum. Connecting the 
conduction and valence bands are the edge states. 



Topological aspects: Chern numbers and edged states 

Topological insulators (TI) are unconventional states of matter [3] that are insu- 
lating in the bulk but conduct along the boundaries. A more general working 
definition of TI, applicable to the case of neutral atoms, is that they are band in- 
sulators with a bulk gap but robust gapless boundary (edge or surface) modes ( 
Fig. [1]). These boundary modes are protected against small deformations in the 
system parameters, as long as generic symmetries such as time-reversal or charge- 
conjugation are preserved and the bulk gap is not closed. Such states arise even in 
systems of non-interacting fermions and are considerably simpler than topological 
states of strongly interacting electrons, such as the fractional QH effect. Despite 
such deceiving simplicity, as we will show in this chapter, the subject is very rich 
and rather intricate, with plenty of open questions. 

We will consider two prime examples of TIs in two dimension, the QH and 
QSH states. A topological insulator is not characterized by any local order pa- 
rameter, but by a topological number that captures the global structures of the 
wave functions. For example, associated with each integer quantum Hall state is 
an integer number, the Chern number . The Chern number also coincides with 
the quantized Hall conductance (in unit of e 2 /h) of the integer QH state. The 
quantization of the Hall conductance is extremely robust, and has been measured 
to one part in 10 9 [ 11 1. Such precision is a manifestation of the topological na- 



ture of quantum Hall states. A quantum Hall state is accompanied by chiral edge 
modes, that propagate along the edge along one direction and their number equals 
the Chern number. Semi-classically, these edge states can be visualized as result- 
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ing from the skipping orbits of electrons bouncing off the edge while undergoing 
cyclotron motion in external magnetic field. The edge states are insensitive to dis- 
order because there are no states available for backscattering, a fact that underlies 
the perfectly quantized electronic transport in the QHE. This well known exam- 
ple demonstrates the "bulk-boundary correspondence", a common theme for all 
topological insulators. 

The QSH state is time-reversal invariant p"3| . It is characterized by a binary 
Z 2 topological invariant v — 1. The edge modes in this case come in the form 
of Kramers pairs, with one set of modes being the time-reversal of the other set. 
In the simplest example, spin up electron and spin down electron travel in op- 
posite direction along the edge. A topologically nontrivial state with v = 1 has 
odd number of Kramers pairs at its edge. Topologically trivial states have u — 0, 
and even number of Kramers pairs, the edge states in this case are not topological 
stable. This is another manifestation of the "bulk-boundary correspondence". Re- 
cent theoretical research has achieved a fairly complete classification, or "periodic 
table", for all possible topological band insulators and topological superconduc- 
tors/superfluids in any spatial dimensions |4]-|6). Several authoritative reviews of 
this rapidly evolving field now exist [[3j|7}[8), which also cover the fascinating sub- 
ject of how to describe TIs using low energy effective field theory with topological 
(e.g., Chern-Simons or 9-) terms. 

In cold atom systems, transport measurement such as Hall conductance is dif- 
ficult. Also, the existence of an overall trap potential instead of a sharp edge for 
a solid state sample, as well as the lack of a precision local probe such as STM, 
makes it hard to probe the edge states. Therefore, identifying topological state 
such as a QH state and extract the Chern number remains a challenging task. Here 
we show that Chern numbers are "hiding" in time of flight images of the cold atom 
systems, a routine measurement in cold atom experiments. 



The Hofstadter model and the Chern number 

For a moderately deep square optical lattice, the dynamics of (spinless) fermionic 
atoms in the presence of an artificial magnetic field is described by the tight bind- 
ing Hamiltonian known as the Hofstadter model fll4| , 

ZJ — ^ 4- p\ p^n + l,m:ri,m ^ _|_ 4- p\ p^n,m+l:n,m _|_ U ^ (\\ 
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Figure 2: The energy spectrum as function of the magnetic flux a per plaquette, 
commonly referred as the butterfly spectrum, for the Hofstadter model with A = 
t y /t x = Irrespectively. 



where c n , m is the fermion annihilation operator at site (n, m), with n/mis the site 
index along the x/y direction, t x and t y are the nearest-neighbor hopping along 
the x and the y-direction respectively. The anisotropy in hopping is described 
by dimensionless parameter A = t y /t x . The phase factor Q^m-.n 1 ,m' is given the 
line integral of the vector potential A along the link from site (n', m') to (n, m), 
(2-Ke/ch) f A • d\. The magnetic flux per plaquete measured in flux quantum is 
labeled a, and represents the strength of magnetic field perpendicular to the 2D 
plane. The energy spectra for A = 1,2 are shown in Fig. [2] At commensurate 
flux, i.e., a is a rational number a = p/q with p and q being integers, the original 
energy band splits into q sub-bands by the presence of magnetic field, and the 
spectrum has q — 1 gaps. For incommensurate flux, namely when a is an irrational 
number, the fractal spectrum contains a hierarchy of sub-bands with infinite fine 
structures. As shown in Fig. [2j varying anisotropy changes the size of the gaps 
without closing any gap. 

In the Landau gauge, the vector potential A x = and A y = —ax and the 
unit cell is of size q and the magnetic Brillouin zone is — n/q < k x < 7i/q and 
— 7r < k y < 7r (the lattice spacing is set to be 1). The eigenstates of the system can 
be written as \l/ nim = e lk ^ n + lk y m ^ !n where tp n satisfies the Harper equation [ 14 1, 



A <+i + e~ lkx V n -i + 2A cos(27ma + k y W n = —Eijj r n . (2) 



e 



Here ip r n+q = ^f n , r = 1, 2, ...q, labels linearly independent solutions. The prop- 
erties of the two-dimensional system can be studied by studying one-dimensional 
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Figure 3: Left: energy spectrum for A = 10, a = 13/21 . The red dots indicate 
dimerized states at the band edges. Middle: the wave function of a state localized 
at a single site, with energy away from the band edge (black), and the dimerized 
wave function that localizes at two distinct sites at a distance of four-lattice spac- 
ing apart for state at band edge near a gap with Chern number 4 (red). Figure on 
the right shows the two point correlation function C(l) that vanishes unless I = 4 
resulting in almost sinusoidal momentum distribution with period 4. 

equation ([2]). For chemical potentials inside each energy gap, the system is in an 
integer quantum Hall state characterized by its Chern number c r and the transverse 
conductivity a xy = c r e 2 / h, where the Chern number of r-fiHed bands is given by, 

c T = Im^ I dk x dk y ^d kx {il) r n yd ky tl) r n . (3) 

r=l ^ n=l 

Here the integration over k x , k y is over the magnetic Brillouin zone. 

Chern-Dimer mapping for anisotropic lattices 

As illustrated in Fig. Q, the hopping anisotropy changes the size of the gaps, 
without ever closing the gap. Hence the topological aspects of the underlying 
states are preserved, as A is continuously tuned. In the limit A 3> 1, the problems 
becomes tractable using degenerate perturbation theory p3|. For A — > oo, the q 
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bands have dispersion E r {k y ) = — 2t y cos(27rar + k y ), r e {l,...,q}. These 
eigenstates are extended along y but localized to single sites as a function of x, 
with wave function i\f n = 8 n>r . These cosine bands shifted from each other inter- 
sect at numerous degeneracy points k*. For finite A, the hopping term proportional 
to t x hybridizes the two states at each crossing, resulting an avoided crossing as 
shown in Fig Q. Near each degeneracy point the system is well described by a 



two-level effective Hamiltonian [ 15 1, 



H(k) = A e [cos(k x £)a x + sm(k x £)a y ] + v e (k y - k*)a z . (4) 

Here a x ,y,z are the Pauli matrices in the space of localized states |r) and |r'), 
and £ = r' — r is the spatial separation between these two localized states. The 



expression for A? and vi are given in Ref. [ 16 1. The eigenstates are equal- weight 
superpositions of \r) and \r'), 

^ = ^(<W + e i/3± M> (5) 

where f3 + = —£(k x + n) + ir and /3_ = — £(k x + ir) are relative phases for the 
upper and lower band edges, respectively. We call such dimerized state "Chern- 
dimers". Starting from Eq. ([5]), it is straightforward to show [15| that the cor- 



responding Chern number is related to the phase factor e lkxe , and simply given 
by c r = I. Thus, for sufficiently anisotropic hopping, dimerized states form at 
the band edges with its spatial extent along x equal to the Chern number of the 
corresponding gap. Such one-to-one correspondence between the Chern number 
and the size of the real space dimer (in the Landau gauge) will be referred to as 
"Chern-dimer mapping". While calculation of Chern number in this limit is not 



new, this mapping has not been elucidated and emphasized until our work [ 16 1. 

The formation of Chern-dimers of size I = c r at the edges of the r-th gap 
implies a "hidden spatial correlation". The correlation function 

C(l) = J2«m^ + l,m) (6) 



will peak at / = £ = c r . For chemical potential in the r-th gap, as A — > oo, ( Fig. 
[4]) C(l) asymptotes to a delta function C(l) — > 5(1 — £), due to the contribution 
of Chern-dimers at the lower edge of the r-th gap. Note that the net contribution 
to C(l) from all other dimerized states associated with gaps fully below the Fermi 
energy is zero, because f3+ — (3- = n for each of these crossings. Remarkably, 
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C{1), or more precisely its Fourier transform n(k x ), is directly accessible in the 
time-of-flight (TOF) image, an routine observable in cold atom experiments, 

In the time-of-flight experiments, the artificial gauge field is first turned off 
abruptly. This introduces an impulse of effective electric field that converts the 
canonical momentum of atoms (when gauge field was on) into the mechanical 
momentum (these subtleties are discussed in [[T6j). Then the optical lattice and 
trap potential are turned off, the density of the atomic cloud n(x) is imaged af- 
ter a period of free expansion £tof- The image is a direct measurement of the 
(mechanical) momentum distribution function with k = mx/htTOF, 

n(k) = ^(ctcv)^-^, (7) 

r,r' 

where the site label r = mx + ny, and for simplicity we neglect the overall 
prefactor and the Wannier function envelope. Thus, from TOF images, the ID 
momentum distribution can be constructed by integrating n(k) over k y , n(k x ) = 
dk y n(k), which is nothing but the Fourier transform of C(l). Thus, in the 
limit of large A, C(l) becomes a delta function, and the ID momentum distribution 
function becomes ^ 

n(k x ) = r ± -cos(c r k x ) + ... (8) 

It assumes a sinusoidal form, with oscillation period given by the Chern number. 



Note the sign depends on which gap the chemical potential lies in [ 1 6 1 , and the off- 
set r comes from the contribution from the filled bands below the Fermi energy. 
In summary, we have arrived at the conclusion that Chern numbers can be simply 
read off from n(k x ) as measured in TOF, in the large anisotropy limit. 



Chern number in time-of-flight images 

We now address the question whether it is still possible to extract c r from TOF 
when away from the large A limit. Note that, while the Chern number does not 
change inside the same gap as A is tuned, n(k x ) is not topologically invariant. 
There is a prior no reason to expect Chern number continues to leave its finger- 
prints in n(k x ). Motivated by these considerations, we have carried out detailed 
numerical studies of momentum distribution in the entire a — A parameter plane, 
going well beyond the "Chern-dimer regime". Fig. ([4]) highlights our results 
with various A but fixed a and Chern number. As A is reduced, n(k x ) gradually 
deviates from the sinusoidal distribution. Yet, the local maxima (peaks) persist. 
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In the limit of isotropic lattice, n(k x ) shows pronounced local peaks. For a wide 
region of parameters, we can summarize our finding into an empirical rule: the 
number of local peaks in n(k x ) within the first Brillouin zone is equal to the Chern 
number. Note that this rule is not always precise, as seen in the Fig. Q for the 
case of Chern number 3 for intermediate values of A where local peaks become a 
broad shoulder-like structure. Despite this, it is valid in extended regions includ- 
ing both large A and A = 1. In fact, given the reasoning above, our empirical rule 
seems to work surprisingly well. But this is an encouraging result for cold atoms 
experiments. From these distinctive features in n(k x ), the Chern number can be 
determined rather easily by counting the peaks, especially for the major gaps with 
low Chern number. 

The Chern-dimer mapping provides a simple, intuitive understanding of n(k x ) 
in the limit of large A. For the isotropic case, there is no dimerization, why are 
there pronounced peaks in n(k x ), and why is the peak number precisely equal to 
the Chern number? To obtain further insight, we considered the low flux limit, 
a — > 0, where the size of the cyclotron orbit is large compared to the lattice 
spacing, and the system approaches the continuum limit. The spectrum then can 
be approximated by discrete Landau levels. For r filled Landau levels, n(k x ) can 
be obtained rather easily, 

n(k x ) = ^(2"i/!)- 1 f^(^)e-« l (9) 

i/=0 

where H v is the Hermite polynomial, and k x is measured in the inverse magnetic 
length ^JeB/hc. n(k x ) is given by the weighted sum of consecutive positive- 
definite function, H^(k x ). The node structure of H u is such that n(k x ) has exactly 
r local maxima (peaks) for r filled Landau levels. This nicely explains the numer- 
ical results for the Hofstadter model at low (but finite) flux a = 1/q, where the 
Chern number c r = r in the r-th gap, and n(k x ) has precisely r peaks (see Fig. [5]). 
Note that this calculation does not directly explain the feature shown in Fig. |4]for 
a = 0.45 and A = 1. As a increases, for [i residing within the same energy gap, 
the distinctive features of n(k x ) are generally retained, so the counting rule is still 
valid. Compared to the low flux regime, however, the local peaks (except the one 
at k x = 0) move to higher momenta. 
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Figure 4: Momentum distribution as anisotropy parameter varies from 1 — 10, 
for a = .45 for chern numbers 1 (top-left), 2 (top-right), 3 (bottom-left), and 4 
(bottom-right). 
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Figure 5: Left: n(k x , /i) in false color for a 50 x 50 lattice with isotropic hopping 
and flux a = 0.05. Right: n(k x ) for chemical potential fi/t=-3A, -2.8, -2.4, -1.8, 
-1.5 respectively (bottom to top). The corresponding Chern number is 1, 2, 3, 4, 
and 5. 

Incommensurate flux: Chern tree 

While it is numerically straightforward to investigate the n{k, x ) or a xy for QHE 
at incommensurate flux a, analytical analysis of the Chern number for this case 
represents a conceptual as well as a technical challenge. In fact, for irrational 
a, the system no longer has translational symmetry, and the notion of magnetic 
Brillouin zone or Bloch wave function lose their meaning. The incommensura- 
bility acts like correlated disorder. The TKNN formula for a periodic system is 
not directly applicable. In addition, the energy spectrum becomes infinitely frag- 
mented (a Cantor set). One expects then a dense distribution of Chern numbers 
even within a small window of energy. Then two questions arise. First, will these 
Chern numbers largely random or are there deterministic patterns governing their 
sequence as ji is continuously varied? Second, what is the experimental outcome 
of n{k x ) or a xy when measured with a finite energy resolution? 

We address some of these questions here for the isotropic lattice (A = 1) us- 
ing the golden-mean flux, a = (\/5 — l)/2, as an example. We will employ a 
well known trick, and approach the irrational limit by a series of rational approx- 
imation a n = F n /F n+ i, where F n is the n-th Fibonacci numbers. Explicitly, the 
series 3/5, 5/8, 8/13, 13/21, 21/34, 34/55, 55/89... approaches the golden-mean 



11 



with increasing accuracy. At each stage of the rational approximation, the energy 
spectrum and the corresponding sequence of Chern numbers can be determined. 
Higher Chern numbers are typically associated to smaller gaps in the asymptoti- 
cally fractal spectrum. At first sight, the Chern number associated with the r-th 
gap seems an irregular function of r as shown in Fig. ^ for a = 34/55. This is 
in sharp contrast with the case of a — 1/q, where c r = r for the r-th gap. 

Now we show analytically that there actually exists a deterministic pattern 
underlying the sequence of Chern numbers. This pattern follows from the so- 
lutions of the Diophantine equation (DE) [[121, r = c r p + s r q, where r is the 
gap index, c r is the Chern number, at each level of the rational approximation, 
p = F n , q = F n+ i. We first classify all the Chern numbers into two classes, the Fi- 
bonacci Chern numbers and the non-Fibonacci Chern numbers. All the solutions 
to the DE, which determine the Chern number, can be exhausted by following two 
simple rules, labeled as I, II and explicitly stated below. 

(I) Fibonacci Chern numbers occur at gaps labeled by a Fibonacci index 
r = F m : for a = F n _ x lF n , c r=Fm = -(-l) n_m F n _ m . Here r = 1 corresponds 
to the first gap near the band edge with Chern number F„_ 2 = (— l) n/ (<? ~ p), 
while the last Fibonacci gap near the band center has a Chern number (— l) n f 
where nf counts the number of Fibonacci numbers less than \q/2\. The Fibonacci 
Chern numbers increases monotonically in magnitude as one moves from the cen- 
ter towards the edge of the band (see Fig. [6]). 

(II) The rest are the non-Fibonacci Chern numbers. The Chern number for the 
r-th gap and its neighboring gaps (r ± 1) obey the recursion relation, 

C r = (C r+ 1 + C r _i)/2, 

c r+ i = c r + (-l) n f(q-p), 
c r _i = c r - (-l) n f(q-p), 

Following these rules, we can construct an exhaustive list of all Chern num- 
bers at each level of rational approximation. We have numerically checked that 
the Chern numbers computed from the TKNN formula are consistent with the pre- 
dictions from the DE above. This procedure yields a hierarchy of Chern numbers, 
which can be neatly organized into a tree structure, thanks to the recursion relation 
above. A subset of this tree hierarchy is shown in Fig. [6} Most notably, we find 
the average of the left and right neighboring Chern number is exactly the Chern 
number in the middle. 

This property helps to understand the measurement outcomes of real experi- 
ments, which inevitably involves average of states with energy very close to each 



12 




Figure 6: The Chern tree at golden-mean flux showing four consecutive rational 
approximations: a = 3/5,8/13,13/21,34/55 starting with Chern number — 1. 
The horizontal axis is not to scale and only shows the correct order of Chern 
numbers. Top-left figure shows Fibonacci (red) and non-Fibonacci (blue) Chern 
numbers vs the gap index r for a = 34/55. Top-right shows momentum distribu- 
tion for three different topological states characterized by Chern numbers 1 (red), 
7 (blue) and —9 (green). Figure illustrates self-averaging aspect of the momentum 
distribution: states 7 and —9 that reside in close proximity to the state with Chern 
number unity have ripples that tend to cancel each other. 
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Figure 7: Left: the quantized conductance c r = a xy /(e 2 /h) versus the chemical 
potential /i for A = 1 and a = (\/5 — l)/2. Right: the blow up for /x close to 
zero showing an extended plateau at c r = 4. c r is computed for a four-terminal 
quantum Hall bar of width W = 20a (inset). The device is connected to four 
semi-infinite leads of the same material with inelastic scattering rate 0.0051 




other. In the context of cold atoms, due to the overall trap potential, the local 
chemical potential slowly varies in space. These local contributions sum up in 
TOF images according to the local density approximation. The self-averaging 
property to some degree ensures that such average yields definite result, and the 
infinitely fine details in the energy spectrum are largely irrelevant. This aspect is 
illustrated in Fig. [6] for the Chern tree shown in Fig. [6} The momentum distribu- 
tion n(k x ) intheaverage sense for a = 3/5,8/13, 13/21, 34/55 will be essentially 
the same. 

Similar consideration applies to a xy . One may see a robust plateau, which 
may not be perfectly quantized, in a region of chemical potential due to the self- 
averaging. We have numerically computed a xy using the Landauer-Buttiker for- 
malism [ 18 1 for a mesoscopic quantum Hall bar (inset of Fig. [7]) at golden-mean 
flux and connected to four semi-infinite leads. The transmission coefficients 
are obtained by solving for the Green's function of the sample and the level width 
functions of each leads. The transverse conductivity at golden-mean flux is shown 
in Fig. |7J One clearly sees a conductance plateau with c r = 4, in addition to 

Gr> — X , 2. . 3 . 
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Quantum spin Hall effect for ultracold atoms 



In solid state, the QSH effect is closely related to spin-orbit coupling. In a recent 
study p7| , it was shown that the effect of spin-orbit coupling can be achieved for 
cold atoms by creating SU{2) artificial gauge fields. 

The theoretical model that captures the essential properties of this system is 
given by a Hamiltonian describing fermions hopping on a square lattice in the 
presence of a SU(2) gauge field 

T~L = t ^ ^ Cm+l.n^ X Cm,n ^m,n+l^ y ^m,n *i" h.C. + A s tag( 1) cJ n n C m n (10) 

mn 

where c m n is the 2-component field operator defined on a lattice site (m, n), and 
Ktag is the onsite staggered potential along the x-axis. The phase factors Q m , n 
corresponding to SU(2) gauge fields are chosen to be, 

9 X = 2vT7cr x , 

6 y = 27txaa z , (11) 

where cr x z are Pauli matrices. The control parameter 7 is a crossover parameter 
between QH and QSH states as 7 = describes uncoupled double QH systems 
while 7 = 7r/2 corresponds to the two maximally coupled QH states. It should be 
noted that the Hamiltonian ( fit)) ) is invariant under time-reversal. The two compo- 
nents of the field operators correspond in general to a pseudo-spin 1/2, but could 
also refer to the spin components of atoms, such as 6 Li, in a F = 1/2 ground 
state hyperfine manifold. For 7 = case, spin is a good quantum number and 
the system decouples into two subsystems which are respectively subjected to the 
magnetic fluxes a and —a. 



For isotropic lattices, this Multiband system has been found to exhibit [ 17 1 a 
series of quantum phase transitions between QSH, normal insulator and metallic 
phases in two-dimensional parameter space for various values of chemical poten- 
tial. Figure Q shows the landscape of the phases near half-filling for a = 1/6. 
The transitions from topologically nontrivial quantum spin Hall (QSH) to topo- 
logically trivial normal band insulator is signaled by the disappearance of gapless 
Kramer-pair of chiral edge modes. 

Motivated by the fact that anisotropic lattices, reveal a novel type of encoding 
of nontrivial topological states in the form of "Chern-Dimers", for U(l) gauge 
systems, we present our preliminary results of QSH states for anisotropic lattices. 
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Figure 8: Shows the bulk (red) and surface spectrum (black) , near half-filling, 
for a = 1/6, 7 = .25, A = 1 , and X stag = .5 ( left) and X stag = 1.5 ( right). 
The presence of gapless edge modes for \ stag = -5 and their absence for \ stag = 
1.5 shows that these two respectively correspond to topologically non-trivial and 
trivial phases. The bottom caption shows the corresponding Phase diagram in 
7 - >^stag plane . 
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Figure 9: Figure illustrates quantum phase transition from topologically non- 
trivial state ( indicated by gapless edge modes ) to normal insulating state, for 
7 = .25, A = 5 , and X stag = 5 ( top-left) and \ sta g = 10 (right). In the bottom 
caption, Green and blue show the local density of the two-component state at all 
sites. For normalized spinor wave function, a density of .25 shows the existence 
of a dimerized state for each component of the spinor. 
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As shown in the Fig. [9j anisotropic lattices subjected to SU (2) gauge fields 
exhibit quantum phase transitions from QSH to normal insulating state by tuning 
the staggered potential. In contrast to the QH system, the avoided crossings in 
QSH system need not be isolated as seen near k y = 1/3. However, there are 
states, such as those near k y = where one sees an isolated avoided crossing. 

Our preliminary investigation as summarized in figure Q shows that for 
anisotropic lattices, in QSH phase, both components of the spinor wave function 
are dimerized near an isolated avoided crossings. It is the existence of this double 
or spinor-dimer that distinguishes QSH state from a normal insulator. In other 
words, each component of the wave function of topologically non-trivial insulat- 
ing phase localizes at two distinct sites, which is in contrast to the topologically 
trivial phase, where the-components of the wave functions localize at a single site. 
For normalized states, the existence of a double-dimer can be spotted by density 
approaching .25 as illustrated in the Fig. [9] The distinction between the QSH and 
normal insulator, characterized by spinor-dimer and its absence suggests that such 
topological phase transitions may be identified by time of flight images. 

Open questions and outlook 

Identifying topological states in time of flight images would be a dream measure- 
ment for cold atom experimentalists. Therefore, our result that the momentum 
distribution carries fingerprints of nontrivial topology is important. 

Our studies raise several open questions regarding topological states and its 
manifestations. Firstly, our detailed investigation suggests that the features of 
topological invariant (Chern number) in momentum distribution are quite robust. 
Is there a rigorous theoretical framework that can establish this connection even 
for parameters beyond Chern-dimer regime? Secondly, the sign of the Chern num- 
ber can not be directly determined from the momentum distribution. It remains an 
open question whether some other observable can detect its sign. 

Another natural question is: how generic is the existence of dimerized states 
in topological states of matter? Somewhat analogous to edge states, the size of 
the dimerized states in our example is a robust (topologically protected) quantity. 
Furthermore, are edge states and dimerized states related to each other? Naively, 
one expects a dimer of size N in bulk will result in iV-possible types of defects 
at open boundaries. These "defects" quite possibly will manifest as edge states 
inside the gap. It is conceivable that the topological aspect of the dimers of size 
N is intimately related to iV-edge states in QH systems. Further work is required 
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to test these speculations, even in special cases. 

To investigate the generality of the dimer framework in topological states, we 
have studied QSH states for anisotropic lattices. Our preliminary studies suggest 
that the topological phase transition from QSH to normal insulating state is ac- 
companied by the absence of dimerized states and hence may be traceable in time 
of flight images. Detailed characterization of Z 2 invariance within the framework 
of spatially extended objects, however remains appealing but open problem. 

The subject of topological states of matter is an active frontier in physics. 
Lessons learned from their realizations in ultracold atoms will help to advance 
the fundamental understanding, and inspire new material design or technological 
applications. Ultracold topological superfluids, Mott insulators with nontrivial 
topology, and non-equilibrium topological states are most likely to bear fruit in 
reasonably near future. 
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